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Abstract 

We discuss direct CP violation in the standard model by giving a new estimate of e'/e in 
kaon decays. Our analysis is based on the evaluation of the hadronic matrix elements of the 
AS = 1 effective quark lagrangian by means of the chiral quark model, with the inclusion of 
meson one- loop renormalization and NLO Wilson coefficients. Our estimate is fully consis- 
tent with the AI —1/2 selection rule in K — ► irir decays which is well reproduced within the 
same framework. By varying all parameters in the allowed ranges and, in particular, taking 
the quark condensate — which is the major source of uncertainty — between (—200 MeV) 3 
and (-280 MeV) 3 we find 

-5.0 x 10~ 3 <e'/e< 1.4 x 10~ 3 . 

Assuming for the quark condensate the improved PCAC result (qq) = —(221 ± 17 MeV) 3 
and fixing Aq^ d to its central value, we find the more restrictive prediction 

e'/s= (4 ±5) x 10~ 4 , 

where the central value is defined as the average over the allowed values of Im At in the 
first and second quadrants. In these estimates the relevant mixing parameter Im At is self- 
consistently obtained from e and we take mj? olc = 180 ± 12 GeV. Our result is, to a very 
good approximation, renormalization-scale and 75-scheme independent. 
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1 Introduction 



The real part of e'/e measures direct CP violation in the decays of a neutral 
kaon in two pions. It is a fundamental quantity which has justly attracted a 
great deal of theoretical as well as experimental work. Its determination would 
answer the question of whether CP violation is present only in the mass matrix 
of neutral kaons (the superweak scenario) or is instead at work also directly in 
the decays. 

On the experimental front, the present results of CERN (NA31) |l| 

Re (e'/e) = (23 ± 7) x 1(T 4 (1.1) 

and Fermilab (E731) @ 

Re (e'/e) = (7.4 ± 6.0) x 1(T 4 (1.2) 

are tantalizing insofar as the superweak scenario cannot be excluded and the 
disagreement between the two outcomes still leaves a large uncertainty. The 
next generation of experiments — presently under way at CERN, Fermilab and 
DA$NE — will improve the sensitiveness to 1 x 10~ 4 and hopefully reach a 
definite result. 

On the theoretical side, much has been accomplished, although the intrinsic 
difficulty of a problem that encompasses scales as different as m t and m n 
weights against any decisive progress in the field. 

A fundamental step was recently covered by the Munich || and Rome 
groups who computed the anomalous dimension matrix of the ten relevant 
operators to the next-to-leading order (NLO) in two 75-schemes of dimensional 
regularization: 't Hooft-Veltman (HV) and Naive Dimensional Regularization 
(NDR). This computation has brought the short- distance part of the effective 
lagrangian under control. 

The residual (and, unfortunately, largest) uncertainty is due to the long- 
distance part of the lagrangian, the computation of which implies evaluating 
the hadronic matrix elements of the quark operators. It is here that the non- 
perturbative regime of QCD is necessarily present and our understanding is 
accordingly blurred. 



1 



At present, there exist two complete estimates of such hadronic matrix 
elements performed by the aforementioned groups, and recently updated in 
ref. for the lattice (at least for some of the operators), where the value 

Re (e'/e) = (3.1 ± 2.5 ± 0.3) x 1(T 4 (1.3) 

is found, and ref. || for the 1/N C approach (for all ten operators) improved by 
fitting the AI = 1/2 rule, where e'/e is estimated to be within the range 

- 2.5 x 1(T 4 < Re (e'/e) < 13.7 x 1(T 4 . (1.4) 

The smaller error in eq. ( |1.3|) originates in the Gaussian treatment of the 
uncertainty in the input parameters with respect to the flat la error included 
in eq. (|L4]). 

Both groups seem to agree on the difficulty of accommodating within the 
standard model a value substantially larger than 1 x 10~ 3 . This unexpectedly 
small value is the result of the cancellation between gluon and electroweak 
penguin operators f7|. If that is actually the case, it is somewhat disappointing 
that the presence of direct CP violation in the standard model turns out to 
be hidden by an accidental cancellation that effectively mimics the superweak 
scenario. 

It seemed to us that a third, independent estimate of e'/e was desirable 
and we have taken the point of view that a reliable evaluation of the hadronic 
matrix elements should first provide a consistent picture of kaon physics, start- 
ing from the CP-conserving amplitudes and, in particular, by reproducing the 
AI =1/2 selection rule, which governs most of these amplitudes as well as the 
quantity e'/e itself. We also felt that the same evaluation should pay partic- 
ular attention to the problem of achieving a satisfactory 75-scheme and scale 
independence in the matching between the matrix elements and the Wilson 
coefficients, the absence of which would undermine any estimate. 

In a preliminary work [JSJ, we studied e'/e within the chiral quark model 
(xQM) in a toy model that included the leading effect of the two most 
important operators, and verified that the 75-scheme independence could be 
achieved. 
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In a recent paper [nj, hereafter referred as I, we have completed the study 



of the hadronic matrix elements of all the ten operator of the AS = 1 effective 



quark langrangian by means of the xQM an d verified in []TT] , hereafter referred 
as II, that the inclusion of non-perturbative 0(a s N c ) corrections and one- loop 
meson renormalization provided an improved scale independence and, more 
importantly, a good fit of the AI =1/2 selection rule. 

These results put us in the position to provide a new estimate of e'/ e that 
is independent of the existing ones and that contains new features that, in our 
judgment, makes it more reliable. 

We summarize here such features. Our estimate takes advantage, as the 
existing ones, of 

• NLO results for the Wilson coefficients; 

• up-to-date analysis of the constraints on the mixing coefficient Im X t . 
Among the new elements introduced, the most relevant are 

• A consistent evaluation of all hadronic matrix elements in the xQM (in- 
cluding non-perturbative gluon condensate effects) in two schemes of 
dimensional regularization; 

• Inclusion in the AS = 1 chiral lagrangian of the complete bosonization 
0{p 2 ) of the electroweak operators Q 7 and Q 8 . Some relevant 0{p 2 ) 
terms have been neglected in all previous estimates; 

• Inclusion of the meson-loop renormalization and scale dependence of the 
matrix elements; 

• Consistency with the AI = 1/2 selection rule in kaon decays; 

• Matching-scale and 75-scheme dependence of the results below the 20% 
level. 

Even though our framework enjoys a high degree of reliability, any estimate 
of e'/e necessarily suffers of a systematic uncertainty that cannot be easily 
reduced further. We find that it is mainly parameterized in terms of the 
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value of the quark condensate, the input parameter that dominates penguin- 
diagram physics. For this reason, we discuss first a inclusive estimate based on 
a conservative range of (qq), as well as the variations of all the other inputs: 
m t , lm\ t (which depends, beside m t and m c , on Bk and other mixing angles) 
and Aqcd- Such a procedure provides us with the range of values for e'/e that 
we consider to be the unbiased theoretical prediction of the standard model. 
Unfortunately, this range turns out to be rather large, spanning, as it can be 
seen in the abstract, from —5 x 10~ 3 to 1.4 x 10~ 3 . On the other hand, it is 
as small as we can get without making some further assumptions on the input 
parameters — assumptions that all the other available estimates must make as 
well. 

In order to provide such a more restrictive estimate, we have chosen the 
improved PCAC prediction for the quark condensate and fixed Aqq D to its 
central value. This reasonable, but nevertheless arbitrary choice allows us to 
give the second, and more predictive estimate reported in the abstract. It is 
the latter that should be compared with the current estimates, while, at the 
same time, bearing in mind also the former unrestricted range as a realistic 
measure of our ignorance. 

Such uncertainty notwithstanding, we agree in the end with the main point 
of ref . H , namely that it is difficult to accommodate within the standard model 
a value of e'/e larger than 1 x 10~ 3 . In fact, if our analysis points toward a 
definite prediction, it points to even smaller values, if not negative ones. This 
can be understood not so much as a peculiar feature of the xQM prediction 
as the neglect in other estimates of a class of contributions in the vacuum 
saturation approximation (VSA) of the matrix elements of the electroweak 
operators. This problem is discussed in detail in I. These new contributions 
are responsible for the onset of the superweak regime for values of m t less than 
200 GeV. In our computation, it is the meson renormalization that in the end 
brings back e'/e around zero or positive values. 

The outline of the paper is the following. In section 2 we write the effec- 
tive quark lagrangian, discuss the short-distance input parameters and give 
the Wilson coefficients. Section 3 contains a brief discussion of the xQM eval- 
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uation of the hadronic matrix elements and their corresponding meson-loop 
renormalization. In section 4 we discuss the values of the input parameters 
and in section 5 the effective factors B^s that give the comparison between the 
VSA and the xQM evaluation of the hadronic matrix elements. We begin in 
section 6 our discussion of e'/e by first studying the 75-scheme independence 
and then the contribution of each operator taken by itself. In section 7, we 
give our estimate as a function of the most important input parameters in a 
series of figures and one table. The numerical value of all input parameters 
are collected in a table in the appendix. 

2 The Quark Effective Lagrangian and the NLO Wilson Coef- 
ficients 

The quark effective lagrangian at a scale \x < m c can be written as |T2| 
G, 



= -^EaWGu). (2.1) 

The Qi are four-quark operators obtained by integrating out in the standard 
model the vector bosons and the heavy quarks t, b and c. A convenient and 
by now standard basis includes the following ten quark operators: 



Qi = 


(S«M/3) V -A C"/3^a) V -A ' 


Q 2 = 




$3,5 — 


(^)v-A^ (^)vtA ' 


$4,6 — 




$7,9 = 


| ( sc 0v-A ^1 ( < m)v±a ' 


$8,10 = 





(2.2) 



where a, [3 denote color indices (a,/3 = 1, . . . ,N C ) and e q are quark charges. 
Color indices for the color singlet operators are omitted. The labels (V ± A) 



5 



refer to 7^(1 ±75). We recall that Qi,2 stand for the W-induced current- 
current operators, Q^-q for the QCD penguin operators and Q7-10 for the 
electroweak penguin (and box) ones. 

The functions Zi{p) and are the Wilson coefficients and Vij the Koba- 
yashi-Maskawa (KM) matrix elements; r = —VtdV t */V u dV* s . Following the 
usual parametrization of the KM matrix, in order to determine s'/e, we only 
need the yi(fi), which control the CP- violating part of the amplitudes. 

The size of the Wilson coefficients at the hadronic scale (~ 1 GeV) depends 
on a s and the threshold masses rriw, rrib and m c . In addition, the penguin 
coefficients yi(^i) depend on the top mass via the initial matching conditions. 

The recent determination of the strong coupling at LEP and SLC gives [ 13| 

a s (m z ) = 0.119 ±0.006, (2.3) 

which corresponds to 

J^qcd = 350 ± 100 MeV . (2.4) 

We will use the range in eq. ( |2.4| j for our numerical estimate of e'/e. 
For m t we take the value |L4j 

m P° le = 180 ± 12 GeV (2.5) 

The relation between the pole mass M and the MS running mass m(/i) is given 
at one loop in QCD by JTT 



m(M) = M(q 2 = M 2 ) (l-t^l^ , (2.6) 

For the running top quark mass, in the range of a s considered, we obtain 

m t (mr' c ) ~ 172 ± 12 GeV (2.7) 

which, using the one-loop running, corresponds to 

m t {m w ) ~ 183 ± 14 GeV , (2.8) 

which is the value to be used as input at the mw scale for the NLO evolution of 
the Wilson coefficients. In eq. Q2.8|) we have averaged over the range of J^qcd 
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given in eq. (2.4). We have explicitly checked that taking \i = mf ole as the 
initial matching scale in place of m\y and using correspondingly fn t (m^ ole ) the 
electroweak Wilson coefficients at [A — 1 GeV remain stable up to the percent 
level, while the variation of the relevant gluon penguin coefficients stays below 
15%. The stability worsen by keeping the top mass fixed, while varying the 
matching scale. 

For rrib we take the value 

m p b ole = 4.8 GeV (2.9) 

which falls in the range 4.5 — 4.9 given in []TJJ, and for m c 

mP° lc = 1.4 GeV (2.10) 



which is in the range 1.2 — 1.9 GeV quoted in ||14|| . These are the quark 
threshold values we use in evolving the Wilson coefficients down to the 1 GeV 
scale. We have checked that varying m^ ole within the 4.5—4.9 GeV range affects 
the final values of the Wilson coefficients at the 0.1% level, while varying the 
charm pole mass between 1.2 and 1.9 GeV affects the results at the 15% level 
at most. 

Even though not all the operators in eq. (|2.2j ) are independent, this basis 
is of particular interest for the present numerical analysis because it is that 
employed for the calculation of the Wilson coefficients to the NLO order in a s 
and aell 

In tables 1 and 2 we give explicitly the Wilson coefficients yi(fi) of the ten 
operators at the scale /i = 1.0 and 0.8 GeV, respectively, in the HV and NDR 
schemes. In the xQM the chiral symmetry breaking scale A x turns out to be 
about 0.8 GeV. This sets a preferential scale for the matching of the hadronic 
matrix elements to the Wilson coefficients. We have checked that the QCD 
perturbative expansion is under control. In fact the difference between LO and 
NLO results for all physical amplitudes considered — both real and imaginary 
parts — at fi = 0.8 GeV remains always smaller than 30%. 

In order to test the fi independence of e'/e we vary the matching scale 
between 0.8 and 1.0 GeV, the highest energy up to which we trust the chiral 
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A (4) 


250 MeV 


350 MeV 


450 MeV 




0.113 


0.119 


0.125 


HV 


2/3 


(0.0007) 0.0301 


(0.0007) 0.0390 


(0.0007) 0.0509 


2/4 


(0.0011) -0.0513 


(0.0012) -0.0610 


(0.0012) -0.0723 


2/5 


(-0.0004) 0.0137 


(-0.0004) 0.0163 


(-0.0004) 0.0209 


2/6 


(0.0011) -0.0766 


(0.0012) -0.103 


(0.0012) -0.144 


Vila 


(0.172) -0.0115 


(0.172) -0.0103 


(0.172) -0.0083 


Vs/a 


(0) 0.167 


(0) 0.230 


(0) 0.328 


y 9 /a 


(-1.19) -1.71 


(-1.19) -1.83 


(-1.19) -2.00 




(0) 0.750 


(0) 0.859 


(0) 1.10 


NDR 


2/3 


(0.0017) 0.0268 


(0.0018) 0.0336 


(0.0018) 0.0416 


2/4 


(-0.0019) -0.0491 


(-0.0021) -0.0574 


(-0.0022) -0.0660 


y5 


(0.0007) 0.0031 


(0.0007) -0.0028 


(0.0007) -0.0165 


ye 


(-0.0019) -0.0849 


(-0.0021) -0.119 


(-0.0022) -0.178 




(0.149) -0.0119 


(0.149) -0.0118 


(0.136) -0.0127 


ys/a 


(0) 0.153 


(0) 0.212 


(0) 0.304 


y 9 /a 


(-1.22) -1.71 


(-1.22) -1.83 


(-1.22) -1.99 


yio/a 


(0) 0.674 


(0) 0.843 


(0) 1.07 



Table 1: NLO Wilson coefficients at fi = 1.0 GeV for m t {m w ) = 183 GeV, 
which corresponds to m? t ole = 180 GeV. The values of the coefficients at /i — 
mw are given in parenthesis (a = 1/128). In addition one has 3/1,2 (/■*) = 0. 
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A (4) 


250 MeV 


350 MeV 


450 MeV 




0.113 


0.119 


0.125 


HV 


V3 


(0.0007) 0.0338 


(0.0007) 0.0456 


(0.0007) 0.0602 




(0.0011) -0.0522 


(0.0012) -0.0626 


(0.0012) -0.0741 


2/5 


(-0.0004) 0.0140 


(-0.0004) 0.0192 


(-0.0004) 0.0397 


ye 


(0.0011) -0.0904 


(0.0012) -0.137 


(0.0012) -0.242 




(0.172) -0.0131 


(0.172) -0.0111 


(0.172) -0.0039 


Vs/a 


(0) 0.209 


(0) 0.327 


(0) 0.583 




(-1.19) -1.82 


(-1.19) -2.04 


(-1.19) -2.45 




(0) 0.835 


(0) 1.14 


(0) 1.66 


NDR 


ys 


(0.0017) 0.0294 


(0.0018) 0.0373 


(0.0018) 0.0422 


2/4 


(-0.0019) -0.0493 


(-0.0021) -0.0569 


(-0.0022) -0.0603 


y5 


(0.0007) -0.0014 


(0.0007) -0.0167 


(0.0007) -0.0708 


ye 


(-0.0019) -0.104 


(-0.0021) -0.171 


(-0.0022) -0.353 


yi/a 


(0.149) -0.0138 


(0.149) -0.0156 


(0.149) -0.0274 


ys/a 


(0) 0.189 


(0) 0.294 


(0) 0.511 


y 9 /a 


(-1.22) -1.81 


(-1.22) -2.03 


(-1.22) -2.42 


yw/a 


(0) 0.819 


(0) 1.11 


(0) 1.60 



Table 2: Same as in Table 1 at fi = 0.8 GeV. 



loop corrections computed in I. We find that, in spite of the fact that some 
of the Wilson coefficients vary in this range by up to 50%, the matching with 
our matrix elements reduces the /z-dependence in e'/e below 20% in most of 
the parameter space. We consider this improved stability a success of the 
approach. 
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3 The Hadronic Matrix Elements 



In paper I we have computed all hadronic matrix elements of the effec- 
tive quark operators in eq. (|2.2|) in the framework of the xQM. The matrix 
elements are obtained by the integration of the constituent quarks by means 
of dimensional regularization. The loop integration leads to results that de- 
pend on the scheme employed to deal with 75 but are scale independent. The 
renormalization-scale dependence is introduced in our approach by the meson- 
loop renormalization of the amplitudes, as explained in I. The meson-loop 
corrections together with the gluon-condensate contributions are the most rel- 
evant ingredients in reproducing the AI = 1/2 selection rule in K — > tttt 
decays, (as discussed in II). 

The xQM results are expressed in a double power expansion on M 2 /A^ 
and p 2 /A 2 , where M is a dimensionful parameter of the model which is not 
determined (generically, it can be interpreted as the constituent quark mass in 
mesons) and p is a typical external momentum. 

The value of M is constrained |15[ by experimental data on the decay of 
7T° and 77 to be 

M = 223±9MeV (3.1) 

(and M = 243 ± 9 MeV if higher order corrections are included). The value of 

M = m p /Vl2 ~ 222 MeV (3.2) 

is found by vector-meson-dominance estimates. Finally, in a recent fit of all 
input parameteres of the extended Nambu-Jona-Lasinio model [16fl, it was 
found a value of 

M~ 200 MeV. (3.3) 

While we could simply take these values and thus make the xQM predictive, 
our approach also allows for a self-consistent determination of a range for M 
that can be compared to the above values. 

The idea is that in physical observables the 75-scheme and /i-dependences 
of the matrix elements should balance the corresponding dependences of the 
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NLO Wilson coefficients . In I we have constructed the complete 0(p 2 ) chiral 
representation of the lagrangian in eq. ( j2.1|) , where the local quark operator 
Qi is represented by a linear combination of bosonic operators B a , namely 
Qi J2 a G a (Qi)B a . The effective quark lagrangian is therefore replaced by 
the following chiral representation 

£f =1 = £ Ci(Ai) G a {Qi) B a . (3.4) 

V 2 i,a 

As mentioned above, and discussed at lenght in I, the chiral coefficients 
G a determined via the xQM approach are 75-scheme dependent. While the 
75-scheme dependence arises in the xQM from the integration of the chiral 
fermions, the explicit /i-dependence is entirely due to the chiral loop renormal- 
ization of the matrix elements: 

(6|£f =1 \a) = -% E C Usb) G a {Q t ) (b\B a \a) (jm ld ) , (3.5) 

where we have labeled by a and b the initial and final bosonic states. We remark 
that in our approach the //-dependence of the chiral loops is not cancelled 
by higher order counterterms, as it is usually required in the strong chiral 
lagrangian. 

The renormalization scale dependence is therefore determined order by 
order in the energy expansion of the chiral lagrangian. In this respect there is 
no direct counterpart to the expansion in strong and electromagnetic couplings 
on which the short- distance analysis is based and, accordingly, we refer to the 
explicit //-dependence in the matrix elements as to the long-distance (LD) or 
"non-perturbative" scale dependence. A purely perturbative renormalization 
scale dependence is introduced in the matrix elements by the NLO running 
of the quark condensate, which we include whenever a comparison between 
values at different scales is required. Otherwise, quark and gluon condensates 
are considered in our approach as phenomenological parameters. 

Our aim is to test whether the estimate of observables is consistently im- 
proved by matching the "long-distance" 75-scheme and \i dependences so ob- 
tained with those present in the short- distance analysis (in particular we iden- 
tify hsd with hld)- Whether and to what extent such an improvement is 
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reproduced for many observables and for a consistent set of parameters, might 
tell us how well low-energy QCD is modelled in the xQM-chiral lagrangian 
approach that we have devised. 

In II, we have shown that minimizing the 75-scheme dependence of the 
physical isospin 1 = and 2 amplitudes determines a range for the parameter 
M between 160 and 220 MeV. In II, it was also found that the /1 depen- 
dence induced by the Wilson coefficients is substantially reduced by that of 
the hadronic matrix elements. 

These issues become crucial for e'/e where the 75-scheme dependence in- 
duced by the Wilson coefficients determines an uncertainty as large as 80% 
when using the 1/N C hadronic matrix elements (see for instance ref. ||) which 
are scheme independent. 

In the following, for the reader's convenience, we report from I the expres- 
sions for the isospin amplitudes for all ten operators in eq. (|2.2| ): 



The corresponding one- loop meson corrections are denoted by cbQ^iQi)- The 
Clebsh-Gordan coefficients for the isospin projections can be found in I. 
For the HV case we obtain: 



(Qi/0,2 



(2tt, J = 0,2 



\Qi\K°). 



(3.6) 




(3.11) 
(3.12) 



(3.10) 



(3.7) 



(3.8) 



(3.9) 



(Qi)o = x + a (g 4 ) 




(3.13) 



(3.14) 



12 



2 (gg> 



(Qrh 
(Qs)o 



n F + -X + a (Q 7 ) 

/-(gg) 2 1 V2(qq) v ^2 



(3.15) 
(3.16) 



-2 



f 3 

■J K 
(gg) 



M/ 2 



F+ 2F X (l + *<oo>)+ao(G8) 



(Q 8 ) 2 

(Q 9 )o 
(Q 9 ) 2 

(Qio)o 

(Qio) 2 



f 3 



-ix 

2 

^x 

2 



1 - 



N 



(1-5 



2N r 



(GG) 



l + — 1-6 



N 



Hgg) 



2 

V2 



1 - — 1-6 



N 



Hgg) 



X 



1 + 



N 



(l - 6{GG) 



+ a 2 (Qg) 
ao(Qio) 
+ a 2 (Qio) 



where 



and 



4 



X' = X [1-6 



M 2 

A 2 " 

x 



+ 3 77% 



M 2 



(3.17) 

+ a 2 (Q 8 ) (3.18) 
(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 
(3.24) 



5{gg) is given by 

N c (a s GG/n) 



6 



(GG) 



16% 2 f 4 



(3.25) 



It is ( |3.25| ) that parameterizes the non-perturbative part of the computation 
by the contribution of the gluon condensate {a s GG/iT), as discussed in I. 
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The renormalization of / is taken into account by replacing / with the 
one-loop parameter fi in the tree-level amplitudes, which amounts to replacing 
1/f 3 with l/fl multiplied by 



+ 3 u - 



1.18. 



In the NDR case we find: 



(Qi)o 

(Ql>2 

(<5 2 >o 

(Q 2 >2 

(Qs)o 
(Q 4 )o 
(Q 5 >o 

(Q 6 )o 
(Q 7 ) 



-X 



-1 + w (l - 5 {GG) ) 



^X 
3 

l -x 

3 



1 + F (l " 3{GG) 



2 - ^ (l - 5(GG>) 



1 + ^ ( X " <W>) 

= ^(x'-(5 (GG) x)+ao(g3) 
= X' + a (Q 4 ) 

2 (gg) y// , n(n s 



+ a (Qi) 
+ a 2 (Qi) 
+ a (Q 2 ) 
+ a 2 (Q 2 ) 



(gg) 



2v/3 (qq? 



X" + a {Q 6 ) 

o M ^ 



2 (gg) 



(gg)Av 



i (gg) 



X" 



(g7)2 = ^irl 1 - 3 ^ 



^2 (gg) v/2 



(g 8 ) = 2x/3 



(qq)' 



-2 



f 3 

(gg) 



, M*fi \ (qq) 
\qq)Al) Mfl 



X" 



y + w x (i + 5 {GG) 



+ ao(Qs) 



(3.26) 

(3.27) 

(3.28) 
(3.29) 

(3.30) 

(3.31) 
(3.32) 
(3.33) 

(3.34) 



(3.35) 



(3.36) 



(3.37) 
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(Qsh 



(qq)' 
P 



, M*ft 



-V2M Y'-^X(1 + 6 {GG) ) +a 2 (Q, 



2N r 



(Q 9 )o 

(Q 9 >2 

(Qio)o 

(Qio} 2 

where 
X" 



X - — (2X - X' - 5 {GG) X 



2 



N c 



Hgg) 



1 

2 
V2 



2X - X' (1-5 



Nr. 



Hgg) 



02(Qg) 



X 



X 



1 + 



iV, 



(i-t 



(GG) 



X 1-9 



+ «o(Qio) 
+ a 2 (Qio) • 



+ 3 (m 2 K -ml) 



(3.38) 
(3.39) 

(3.40) 
(3.41) 

(3.42) 
(3.43) 



{Qi)2 = for i = 3, 4, 5, 6 in both schemes. 

Of particular interest are the matrix elements (Qe)o and (Qs)2 which dom- 
inate any estimate of e' /e; their leading effect was included in the toy model 
of ref. §. 

The most striking feature concerning the gluon-penguin operators is the 
linear dependence on the quark condensate that is found in the xQM in con- 
trast to the quadratic one of the VSA. This difference explains the different 
weight that these operators have in the two models. 

Concerning the electroweak-penguin operators, as discussed in I, the terms 
proportional to Y — so far neglected in all estimates — give an important contri- 
bution that makes the electroweak-penguin operators larger and, accordingly, 
the cancellation between electroweak and gluon penguins effective even for the 
present values of m*. We shall come back to this point in section 5. 



4 The Input Parameters 

The quark and the gluon condensates are two input parameters of our 
computation. As discussed in I, their phenomenological determination is a 
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complicated question (they parameterize the genuine non-perturbative part of 
the computation) and the literature offers different estimates. 

For guidance, we identify the condensates entering our computation with 
those obtained by fitting the experimental data by means of the QCD sum 
rules (QCD-SR) or lattice computations. 

A review of recent determinations of these parameters, together with a 
justification of the range below, is given in I. Here we only report the ranges 
that we will explore in our numerical analysis. 

For the gluon condensate, we take the scale independent range 

(^GG) = (376 ± 47 MeV) 4 , (4.1) 

7T 

which encompasses the results of recent QCD-SR analysis |18| . While this is 
a crucial input parameter in the physics of the AI = 1/2 rule (see II), it plays 
only a minor role in a penguin-dominated quantity like e'/e. 
For the quark condensate, we consider the range 

- (200 MeV) 3 < (qq) < -(280 MeV) 3 (4.2) 



which includes the central values and the errors of the QCD-SR [19] and lattice 
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The rather conservative range ( |4.2j ) is the one advocated in I. In II it is 
shown that the AI = 1/2 selection rule seems to prefer the upper half of this 
range. 

For comparison, in most estimates of e'/e the PCAC value 

{QQ) W = - _ / x , _ 7 x , (4-3) 

is taken (with 8k equal to zero) and the error range is that of the determination 
of m. 



m s (l GeV) = 178 ± 18 MeV . (4.4) 
This choice gives a quark condensate of 

(qq) = -(261 ±9 MeV) 3 (4.5) 
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at /i = 1.0 GeV that corresponds, via NLO renormalization at fixed Aqq D = 
350 MeV, to 

= -(244±9MeV) 3 (4.6) 



at the matching scale of fi = 0.8 GeV. The quark condensate in ( |4.6|) has an 
uncertainty of only 9 MeV, that is perhaps too small a range to account for 
the actual uncertainty. For instance, if the QCD-SR improved estimate 



8 K = 0.34j£t? (4.7) 

is taken into account, the range (at /i = 1 GeV) to be explored becomes 
comparable to that of 



- (183 MeV) 3 < (qq) < -(256 MeV) 3 , (4.8) 



with a central value much smaller than in 

The range ( |4.8|) is consistent with what one finds by the same PCAC value 
as in ( [4.3|) but with f K and m K replaced by, respectively, f w and m^, and 
m s (fj.) + m d (fj) by 

fn u (l GeV) + fn d {l GeV) = 12 ± 2.5 MeV , (4.9) 

as given in [|23|| . By means of the latter, taken at \x = 0.8 GeV, we find 

(qq) = (-221 ± 17 MeV) 3 . (4.10) 

The range ( |4. 1C| ) suffers of a larger error with respect to that of ( |4.6j ), which 
however does not take into account ([4.7| ) and, accordingly, the much broader 
range ( [4.8|) which is more realistic. 

Because of such uncertainties, we will consider in section 7 two possible 
ranges for the quark condensate: the range (|4.2| ) for our most conservative 
estimate and the improved PCAC result (|4.10| ) for a second, more restrictive 
one. As discussed in the introduction, these two ranges complete each other 
in providing, at the same time, a definite prediction and a gauge of the overall 
uncertainty of the prediction itself. 

In section 7, in order to make the comparison to other estimates easier, we 
will also give our result for the range ( [4.6|) . 
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5 The B { Factors 



Let us introduce the effective factors 

R (0,2) _ \V»/0,2 ( r ,s 

* " (Q;)oJ A ' 1 j 

which give the ratios between our hadronic matrix elements and those of the 
VSA. They £1X6 cL useful way of comparing different evaluations. 

In table 3, we collect the Bi factors for the ten operators. The values of 
the Bi depend on the scale at which the matrix elements are evaluated, the 
input parameters and M; moreover, in the xQM they depend on the 75-scheme 
employed. We have given in table 3 a representative example of their values 
and variations. 

The values of B± and B^ show that the corresponding hadronic ma- 
trix elements in the xQM are, once non-factorizable contributions and meson 

renormalization have been included, respectively about ten and three times 

(2) (2) 

larger than their VSA values. At the same time, B\ and B\ turn out to be 
at most half of what found in the VSA (for the starred entries see the comment 
at the end of the section). These features make it possible for the selection 
rule to be reproduced in the xQM, as extensively discussed in II. 

For comparison, in the 1/N C approach of ref. ||24j| , the inclusion of meson- 
loop renormalization through a cutoff regularization, leads, at the scale of 1 
GeV, to B{ 0) = 5.2, Bf ] = 2.2 and b{ 2) = B% ] = 0.55, a result that is not 

sufficient to reproduce the AI = 1/2 rule. The similarity of the HV values 

(2) (2) 

B\ = B\ = 0.55 obtained in the xQM with the corresponding 1/N C results is 
remarkable, and yet a numerical coincidence, since the suppression originates 
from gluon condensate corrections in the xQM, whereas it is the effect of the 
meson loop renormalization (regularized via explicit cut-off) in the analysis of 



ref. 121 



The values of the penguin matrix elements (Q3) and (Q4) in the xQM 
lead to rather large Bi factors. In the case of Q3, the xQM result has the 
opposite sign of the VSA result and B 3 is negative. This is the effect of the 
large non-perturbative gluon correction. 
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HV 


NDR 


fi = 0.8 GeV 


= 1.0 GeV 


H = 0.8 GeV 


H = 1.0 GeV 


b[ 0} 


10.6 


11.1 


10.6 


11.1 


o(0) 


2.8 


3.0 


2.8 


3.0 




0.52 


0.55 


0.52 


0.55 




0.52 


0.55 


0.52 


0.55 




-2.9 


-3.0 


-3.7 


-3.9 


B 4 


1.8 


1.9 


1.0 


1.1 


B5 = Bq 


1.7 -=-0.61 


1.8 -=-0.64 


1.0 -=-0.38 


1.1 -=-0.40 


Bf ] 


3.0 -=-2.2 


3.3 -=-2.4 


2.9 4-2.2 


3.242.3 


B 8 


3.3 -=-2.2 


3.6 -=-2.4 


3.2 4 2.2 


3.5 -=-2.4 


(0) 


3.9 


4.0 


3.5 


3.6 


B (o) 


4.4 


4.7 


5.6 


5.9 




2.7 -f- 1.5 


3.0 4- 1.5 


2.7 -=-1.4 


2.9 -=-1.5 


n 8 


2.1 -=- 1.4 


2.3 -=- 1.5 


2.1 4 1.4 


2.3 4- 1.5 


o(2) 


0.52 


0.55 


0.52 


0.55 




0.52 


0.55 


0.52 


0.55 



Table 3: The Bi factors in the xQM (including meson-loop renormalizations) 
at two different scales: fi = 0.8 and 1.0 GeV and in the two 75-schemes. We 
have taken the gluon condensate at the central value of eq. ( |4.1|) , while the 
range given for -65-8 corresponds to varying the quark condensate according 
to eq. ( |4.2[ ). The results shown are given for M = 220 MeV. 

Regarding the gluon penguin operator Q 6 (and Q 5 ), we find that the xQM 
gives a result consistent with the VSA (and the 1/N C approach), B 6 (B 5 ) being 
approximately equal to two for small values of the quark condensate and one- 
half at larger values. It is the quadratic dependence (to be contrasted to the 
linear dependence in the xQM) of the VSA matrix element for the penguin 
operators that it responsible for the different weight of these operators at 
different values of the quark condensate. The lattice estimate at \x = 2 GeV 
for these operators gives B5 = Bq = 1.0 ± 0.2 ||. A direct comparison in this 
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case is not possible. 

The electroweak Bi factors are all larger in the xQM than in the VSA, 

(2) 

except for Bg ( that are about 1/2 in the HV and about 0.4 in the NDR 
scheme. For comparison, the lattice estimate at fi = 2 GeV in this case yields 
Bfl = 1.0 ± 0.2 and Bf ] = 0.62 ± 0.10 g. 

The most relevant result for e'/e is the value of B^ which ranges from 1.5 
to 2 times that of the VSA. This increase is due to two independent reasons. 
On the one hand, we found two new terms in the chiral lagrangian that have 
not been included so far in the VSA estimate of the matrix elements. The chiral 
coefficients of these terms are computed in the xQM approach — as discussed 
in detail in I. It is an open question how they can be determined in the VSA 
framework. 

From this point of view, what we have referred to as VSA — and used in 
table 3 as normalization for the Q-j^ operators — is not the complete VSA result. 
The inclusion of the new terms amounts up to a 60% increase of for small 
values of (qq) in the chosen range and down to about 10% for large values; 
smaller effects are found in the case of -B^g. On the other hand, the meson- 
loop renormalization associated with the new chiral terms is large (see I) and 
adds up to reproduce the results shown in table 3. The increase in importance 
of the operator Q$ with respect to Qq turns into a more effective cancellation 
between the two operators for large values of the quark condensate while at 
smaller values the gluon penguin contribution prevails. 

The relations Bf ] = B%\ B^ = b{ 2) and b[ 2) = B% ] hold true in both 75 - 
schemes. These relationships are a reminiscent of those among the operators — 
which are preserved by gluon corrections and meson renormalization. 
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6 Studying e'/e in the xQM 



The quantity e'/e can be written as 



e 2 e Re A 



G f uj 

Im \ t 



n - - n 2 

00 



(6.1) 



where, referring to the AS = 1 quark lagrangian of eq. (|2.1|) , 

n = !><Qi)o (6.2) 

n 2 = J2Vi(Qi)2+U J2yi(Qi)0 Q V+V' ' ( 6 - 3 ) 

i j 

and 

ImA t = Imy t ^;. (6.4) 

The quantity fi^^ includes the effect of the isospin-breaking mixing be- 
tween 7T° and the etas. 

Since Im A„ = according to the standard conventions, the short- distance 
component of e'/e is determined by the Wilson coefficients yj. Following the 
approach of ref. ||, yi(fi) = y% (//) =0. As a consequence, the matrix elements 
of Q\.2 do not directly enter the determination of e'/e. On the other hand, 
in the HV scheme the matrix elements of Q4 Qg and Qiq can be expressed 
in terms of those of Qi^ and Q%. The work of ref. has taken advantage of 
this fact to determine some of the penguin matrix elements, after imposing the 
AI = 1/2 rule. The xQM determination of (Qi,2,3) gives for {Qa) q a result 
that differs substantially from that used in ref. , as we discuss at the end of 
section 6.3. 

We take, as input values for the relevant quantities, the central values given 
in appendix. We thus have 

( ' hxC 349 GeV" 3 , u = 1/22.2, 0^ = 0.25. (6.5) 



2 lei Re A 







The large value in eq. ( |6.5|) for 1/lu comes from the AI = 1/2 selection 
rule. In II we have shown that such a rule is well reproduced by the xQM 
evaluation of the hadronic matrix elements. As the precise values of Re A and 
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oj depend on the choice of the input parameters and of M — the selection rule 
being satisfied within a 20% approximation — we have taken the corresponding 
experimental values. Similarly, the value taken for e is the experimental one. 



6.1 The Mixing Parameter ImA* 



A range for Im At is determined from the experimental value of c as a 
function of m t and the other relevant parameters involved in the theoretical 
estimate. We will use the most recent NLO results for the QCD correction 
factors 771,2,3 which are given in the NDR scheme [25| and vary the AS = 2 
hadronic parameter B^ around the central value obtained in the xQM using 
the same regularization scheme. 

In order to restrict the allowed values of Im \ t we have solved the two 
equations 



£th{B K , \V cb \, \V US \, \V ub \/\V cb \, A QC D,m t ,m c ,r], p) 



2 i 2 
V +P 



1 


v ub 


2 


IK, 


2 


v cb 


2 



(6.6) 
(6.7) 



to find the allowed values of 7/ and p, given m t , m c and [I4"| 



IK 



IK 



ub\ 



us I 

v cb \ 
v cb \ 



(6.8) 
(6.9) 
(6.10) 
(6.11) 



: (2.266 ± 0.023) x 10" 3 
: 0.2205 ±0.0018 
0.041 ±0.003 
0.08 ±0.02. 

For the renormalization group invariant parameter we take the rather 
conservative range 

B K = 0.55 ± 0.25 (6.12) 

that encompasses both the xQM model prediction [^] and other current de- 
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For the NLO order ^-parameters for Aq^ d = 350 MeV and mf ULK '' = 180 
GeV, at /i = m c we find 



(pole) 



77! = 1.36 772 = 0.513 7/3 = 0.446 



(6.13) 



22 



We do not include bounds provided by the quantity Xd of B-physics that 
we find to have a marginal impact in the determination of e' / e once the large 
error in (|6.12|) is taken into account. 

This procedure gives two possible ranges for ImA t ~ ?y|V^ s | | V c b\ 2 , which 
correspond to having the KM phase in the I or II quadrant (p positive or 
negative, respectively). For example, for mf ole = 180 GeV (m t (mw) — 183 
GeV) and A^ D = 350 MeV we find 

1.1 x 10~ 4 < ImAt < 1.9 x 10~ 4 (6.14) 

in the first quadrant and 

0.75 x HT 4 < ImA 4 < 1.9 x 10 -4 (6.15) 



in the second quadrant. For the range of Bk given in eq. (|6.12j) varying all the 



other parameters (including m t and Aqcd) affects the above limits on lm.X t 
by less than 20%. In particular, the upper bound on ImA 4 is stable and it is 
directly related to the maximum value of rj obtained from eq. ( |6.7|) (p = 0). 
The upper bound on Im X t becomes a sensitive function of the input parameters 



only if we consider Bk > 0.5. In other words, we agree with ref. |27j that it 
is the theoretical uncertainty on the hadronic AS = 2 matrix element that 
controls the uncertainty on the determination of ImA t . 

We have included the bounds provided by eqs. (|6.6|) - (|S.7|) in all the follow- 
ing tables and figures. 

6.2 75-scHEME Independence 

In order to fix M, we compare the computation in the HV 75-scheme with 
that in the NDR. Figs. 1, 2 and 3 show how the intersection between the two 
results remains stable as we change the value of the quark condensate. Fig. 4 
shows the change in stability that occurs as we change Aq<I d . 

We find that the values at which 75-scheme independence is achieved 

M ~ 215 — 220 MeV (6.16) 

are quite stable with respect to different values of the matching, the quark and 
the gluon condensates and m t . Smaller values of M are selected for smaller 
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-0.5 



-1 ■ 

~ L5 0.19 0.2 0.21 0.22 0.23 0.24 0.25 

M (GeV) 

Figure 1: 75-scheme dependence of e'/e. The black (gray) line represents the 
HV (NDR) result. The results are shown for (qq) (0.8 GeV) = (-200 MeV) 3 , 
m P° lc = 180 GeV, ImAt = 1.3 x 10~ 4 and A^ D = 350 MeV. The stability is 
obtained at about M = 215 MeV. e'/e is given in units of 10~ 3 . 

values of A QCD (and a correspondingly higher value of e'/e). These results 
are consistent with those found in II for the AI = 1/2 selection rule, where 
stability is achieved in the range M = 160 — 220. They are also consistent 
with the independent estimates discussed in section 3. 

As it is apparent from the figures, the final value of e'/e strongly depends 
on the value of M we take. It is only through the device of requiring 75-scheme 
independence that we are able to reach a definite prediction. This procedure 
has the precious pay-off of providing us with an improved estimate that does 
not suffer of the uncertainty due to the 75-scheme dependence of the NLO 
Wilson coefficients, which may be as large as 80%. 

Figs. 4 and 5 show how the intersection depends on Aq^ d . 
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0.19 0.2 0.21 0.22 0.23 0.24 0.25 
M (GeV) 

Figure 2: Same as Fig. 1 for (qq) (0.8 GeV) = (-240 MeV) 3 
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6.3 Anatomy of e'/e 

It is useful to consider the individual contribution to e'/e of each of the 
quark operators. We have depicted them as histograms, where the grey (black) 
one stands for the contribution before (after) meson-loop renormalization. 
Henceforth all results are given for M = 220 MeV in the HV scheme. 

It is clear from the histograms of Fig. 6, 7 and 8 that the two dominating 
operators are Qq and Q$. Yet, since they give contributions approximately of 
the same size and opposite in sign, the final value turns out to be relatively 
small and of size comparable to that of most of the other operators. This result 
is at the origin the large theoretical uncertainty as well as the unexpected 
smallness of e'/e. 

e'/e 




Q3 Q4 Q5 Q6 Q7 Q8 Q9 Q10 All 



Figure 6: Histograms of the partial contributions to e'/e of the height relevant 
operators for (qq) (0.8 GeV) = (-200 MeV) 3 , mf^ = 180 GeV, ImA t = 1.3 x 
10- 4 and Aqq D = 350 MeV. Gray (black) histo grams represent the contribution 
of each operator without (with) meson-loop renormalization. The last two 
histograms correspond to the sum of all contributions. 

The same histograms serve the purpose of showing that the meson-loop 
renormalizations are crucial not only in the overall size of each contribution but 
also in determining the sign of the final result (see Fig. 7). These corrections 
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e7e 




Q3 Q4 Q5 Q6 Q7 Q8 Q9 Q10 All 

Figure 7: Same as Fig. 6 for (qq) (0.8 GeV) = (-240 MeV) 3 . 



e7e 




Q3 Q4 Q5 Q6 Q7 Q8 Q9 Q10 All 

Figure 8: Same as Fig. 6 for (qq) (0.8 GeV) = (-280 MeV) 3 . 
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are here consistently included in the estimate for the first time. 

The role of the operator Q4 turns out to be marginal in our approach. In 
comparing this result with that of the 1/N C framework || (see also the final 
tables in ref. where we reproduce the individual 1/N C contributions for 
the standard ten operators), it should be recalled that in the above analysis 
the Qi operator is written in terms of Qi, Q 2 and Q 3 and that its values is 
therefore influenced by the Bi factors assigned to the former matrix elements. 
In particular, while B\ and B2 are in ref. || requested to be large in order to 
account for the AI = 1/2 rule, B 3 is assigned the value of 1. Such a procedure 
produces a rather large value for the matrix element of Q4. In our approach, 
we see that in fact also B 3 is large (and negative!) and that Q 4 , once written 
in terms of the other operators, is small, as found in the direct estimate. 
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7 Estimating e'/e 



The preliminary work of the previous sections allows us to estimate e'/e. 
The two most important sources of uncertainty are the quark condensate and 
the value of Im X t . Accordingly, we plot the values of e'/e as a function of 
these two quantities. Fig. 9 and 10 show our estimates, for m t fixed at its 
central value, in the first and second quadrant respectively. As it can be seen 
by inspecting these figures, the larger the value of the quark condensate, the 
swifter is the change in e'/e. 




Figure 9: e'/e in units of 10~ 3 for mf olc = 180 GeV and A { q ] cd = 350 MeV 
as a function of A = lm\ t and the quark condensate (qq) = (qq) 1 ^ in units 
of GeV. Im X t is taken in the first quadrant. Black (grey) squares represent 
positive (negative) values. 

To have an idea of the effect of varying m t , the third major source of 
uncertainty in the input parameters, we have included Fig. 11 where the top 
mass is varied in the given range. 

Fig. 12 shows the stability of our prediction for different matching scales 
/i = 0.8 and 1 GeV (the perturbative running of (qq) is included by taking the 
value of the condensate at \i = 0.8 GeV as the input value and than running 
it to /j, — 1 GeV). The matching-scale dependence is below 20% in most of the 
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0.5 



e'/e 



o 



-0.5 



-1 



0.000075 



0.0001 



0.000125 



X 0.00015 



0.000175 



-0.28 



Figure 10: Same as Fig. 9 with Im X t in quadrant II. 



range, becoming almost 30% only for very large values of the quark condensate. 

In order to provide the reader with a more analytical view, we have also 
collected in Table 4 the numerical results at the varying of all the relevant 
parameters. 

After twelve figures and four tables, we hope to have convinced the reader 
that the quantity e'/e is difficult to estimate with great precision. We think 
that only the order of magnitude can be predicted in a completely reliable 
manner. The reason is very simple: the final value is the result of the cancel- 
lation between two, approximately equal in size, contributions. Accordingly, 
even a small uncertainty will be amplified and we are unfortunately dealing 
with rather large ones. And yet, the shear importance of this quantity impels 
us to provide the best estimate we can. 

By varying all parameters in the allowed ranges and, in particular, tak- 
ing the quark condensate — which is the major source of uncertainty — between 



(-200 MeV) 3 and (-280 MeV) 3 we find 
-27x 10~ 4 <e'/e< 9x 10~ 4 , 



(7.1) 



where we have kept Aq^ d fixed at its central value. A larger range, 
-50 x 10~ 4 < e'/e < 14 x 10" 4 , 



(7.2) 
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Figure 12: Same as Fig. 9 for different matching scales: /x = 0.8 GeV (upper 
surface) and 1.0 GeV (lower surface). 
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A QCD = 250 Mev 


(qq) {1/3) (McV) 


m pol ° (GeV) 


quadrant I 


quadrant II 


-200 


168 


3.7-^6.5 


2.7- 


-6.3 


180 


2.9^-5.3 


2.0- 


-4.9 


192 


2.1 -e- 3.9 


1.4- 


-4.1 


-240 


168 


3.5-^6.3 


2.6- 


-6.1 


180 


1.6 -e- 3.0 


1.2- 


-2.8 


192 


-0.2-7- -0.4 


-0.1 - 


-0.4 


-280 


168 


-0.1 -7- -0.2 


-0.1 - 


- -0.2 


180 


-4.1 -e- -7.5 


-2.9- 


-7.0 


192 


-8.0-7- -15 


-0.5- 


r -16 


A QCD = 350 Mev 


(qq) {1/3) (MeV) 


m pol ° (GeV) 


quadrant I 


quadrant II 


-200 


168 


5.0-^8.5 


3.7- 


-8.7 


180 


3.8^6.9 


2.7- 


-6.6 


192 


2.6-^4.7 


1.7- 


-4.9 


-240 


168 


4.3-^7.4 


3.2- 


-7.6 


180 


1.5-7-2.7 


1.1 - 


-2.6 


192 


-1.3 -e- -2.4 


-0.9- 


-2.5 


-280 


168 


-2.3 -e- -4.0 


-1.7- 


-4.1 


180 


-8.2-7- -15 


-5.7- 


r -14 


192 


-14-7- -26 


-9.4- 


r -27 


Aqq D = 450 Mev 


(qq) {1/3) (McV) 


m pol ° (GeV) 


quadrant I 


quadrant II 


-200 


168 


8.4-7-14 


6.2-7-14 


180 


6.2-7-11 


4.4-e-ll 


192 


4.0-^7.4 


2.7-^6.9 


-240 


168 


6.5-=- 11 


4.8-7-11 


180 


1.5-7-2.7 


l.l-e-2.6 


192 


-3.4-e- -6.2 


-2.3 -e- -5.8 


-280 


168 


-6.9-e- -12 


-5.1 -e- -11 


180 


-17-7- -30 


-12-7- -30 


192 


-27-7- -50 


-18-7- -47 



Table 4: Table of e'/e in units of 1(T 4 . Matching scale // = 0.8 GeV. The 
two values corresponds to, respectively, the lower and upper bounds of Im X t , 
which are determined by consistency with e. 
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is obtained by varying Aqc D as well. 

It should be stressed that the large range of negative values that we obtain 
is a consequence of two characteristic features of our matrix elements: i) the 
enhancement of the size of the electroweak matrix elements (Qsj) due to the 
coherent effects of the additional 0{p 2 ) contributions so far neglected (see dis- 
cussion in sect. 5) and the chiral loop corrections; ii) the linear dependence on 
(qq) of the leading gluon penguin matrix elements compared to the quadratic 
dependence of the leading terms in the electroweak matrix elements, which 
makes the latter prevail for large values of the quark condensate. The effect 
of i) represents an enhancement of the leading electroweak matrix elements by 
a factor two with respect to the vacuum insertion approximation and present 
1/N C estimates (see table 3), while feature ii) is absent in the 1/N C approach, 
the quark condensate dependence being always quadratic. 

To provide a somewhat more restrictive estimate we may assume for the 
quark condensate the improved PCAC result, namely (qq) = —(221 ±17 MeV) 3 
at our matching scale \i = 0.8 GeV, and thus find 



, i 4.5 tf'A x 10 4 quadrant I 
' 3.9 ±i£ x 10" 4 quadrant II. 



The value of e'/e = (4 ± 5) x 10 4 quoted in the abstract is obtained by 
averaging over the two quadrants in eq. (|7.3j). 



The range ft4.10|) for the quark condensate, on which the above estimate 
is based, is not the favorite one by our analysis of the AI = 1/2 selection 
rule in the xQM. The upper half of the more conservative range (|4.2j) seems to 
accommodate more naturally the rule, at least for a constituent mass M ~ 220 
MeV — the value we find by requiring 75-scheme independence of e'/e. For 
large values of the quark condensate the central values of e' / e shift toward the 
superweak regime, and the role of meson loop corrections becomes crucial. 

By taking the quark condensate in the range ( |4.6| ), the QCD-SR improved 
PCAC result, we find 

f 1.4±1;| x 10- 4 quadrant I 
" '' ' ] 1-2 ±2:g x 10" 4 quadrant II. 1 ' ' 
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Actually, for such a range of (qq), negative central values of e'/e in both quad- 
rants are obtained due to the extra terms of the bosonization of the electroweak 
operators Q7 and Q 8 neglected in the previous estimates. Only after the inclu- 
sion of the meson-loop renormalization e'/e turns to the positive central values 
of eq. (|7.4j) . In Fig. 13 we have summarized the present status of the theoretical 
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Theoretical Predictions 



Figure 13: Present status of theoretical predictions and experimental values for 
e'/e (in units of 10~ 4 ). The most recent 1/N C [6] and lattice [5] estimates are 
compared to (a) our unbiased estimate (7.2), (b) our more restrictive estimate 
(7.3). 

predictions for e'/e, compared to the present 1 a experimental results. 



8 Outlook 



Our phenomenological analysis, based on the simplest implementation of 
the xQM and chiral lagrangian methods, takes advantage of the observation 
that the AI = 1/2 selection rule in kaon decays is well reproduced in terms of 
three basic parameters (the constituent quark mass M and the quark and gluon 
condensates) in terms of which all hadronic matrix elements of the AS = 1 
lagrangian can be expressed. 



35 



We have used the best fit of the selection rule to constrain the allowed 
ranges of M, (qq) and (GG) and we have fed them in the analysis of e'/e. 

Nonetheless, the error bars on the prediction of e'/e remain large. This 
is due to two conspiring features: 1) the destructive interference between the 
large hadronic matrix elements of Qq and Q$ which enhances up to an order 
of magnitude any related uncertainty in the final prediction (this feature is 
general and does not depend on the specific approach); 2) the fact that large 
quark-condensate values are preferred in fitting the isospin zero K° — > im 
amplitude at 0(p 2 ) (which is a model dependent result). 

Whereas little can be done concerning point 1) which makes difficult any 
theoretical attempt to predict e'/e with a precision better than a factor two, 
an improvement on 2) can be pursued within the present approach. 

Two lines of research are in progress. On the one hand, we are extending 
the analysis to 0(p 4 ) in the chiral expansion to gain better precision on the 
hadronic matrix elements and to determine in a self-consistent way the poli- 
nomial contributions from the chiral loops; preliminary results indicate that 
the AI = 1/2 rule is reproduced for smaller values of the gluon and quark 
condensates, thus reducing our error bar, in the direction shown by our more 
restrictive estimate. On the other hand, we are studying the AS = 2 sector to 
determine at the same order of accuracy Bk and the K^-Kg mass difference 
by including in the latter the interference with long-distance contributions that 
can be self-consistently computed in the present approach. 

Whether this program is successfull may better determine how much of 
the long range dynamics of QCD is embedded in the present approach and 
increase our confidence on the predictions of unknown observables. 
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Input Parameters 



parameter 


value 


Vud 


0.9753 


V us 


0.2205 ±0.0018 


2 n 

sm Ow 


0.2247 


m z 


91.187 GeV 


my/ 


80.22 GeV 


m b 


4.8 GeV 


m c 


1.4 GeV 


\e\ 


(2.266 ± 0.023) x 10" 3 


\v cb \ 


0.041 ±0.003 


\V u b/V c b\ 


0.08 ±0.02 




180 ± 12 GeV 


B K 


0.55 ± 0.25 


U = /tt+ 


92.4 MeV 


Jk = fx+ 


113 MeV 


m-x — (m v + + m 7r o)/2 


138 MeV 


rriK — mjfo 


498 MeV 


m v 


548 MeV 




2V2^/. 


^77+77' 


0.25 ±0.05 


A (4) 


350 ± 100 MeV 


m u + m d (1 GeV) 


12 ± 2.5 MeV 


m s (1 GeV) 


178 ± 18 MeV 


m 


-(200-^ 280 MeV) 3 


{a s GG/ir) 


(376 ± 47 MeV) 4 



Table 5: Table of the numerical values of the input parameters. 
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